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Abstract. General analytic energy bounds are derived for A-boson systems 
governed by semirelativistic Hamiltonians of the form 

JV JV 
i — l 1 — i<j 

where V(r) is a static attractive pair potential. A translation-invariant model 
Hamiltonian H c is constructed. We conjecture that (H) > (H c ) generally, and 
we prove this for N = 3, and for N = 4 when m = 0. The conjecture is also valid 
generally for the harmonic oscillator and in the nonrclativistic large-m limit. This 
formulation allows reductions to scaled 3- or 4-body problems, whose spectral 
bottoms provide energy lower bounds. The example of the ultrarelativistic linear 
potential is studied in detail and explicit upper- and lower-bound formulas arc 
derived and compared with earlier bounds. 
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1. Introduction 



One-body Hamiltonians H composed of the relativistic expression \/p 2 + m 2 for the 
kinetic energy of particles of mass m and momentum p and of a coordinate-dependent 
static interaction potential V(r), defined as operator sum 

H=y/p> + m? + V(r), 

provide a simple but very efficient tool for the description of relativistically moving 
particles [1, 2, 3]. They have been used, for instance, for the description of hadrons 
as bound states of quarks [4] . One of the advantages of this kind of semirelativistic 
treatment is that its generalization to the many-body problem is straightforward to 
formulate [5]. A semirelativistic Hamiltonian for a system of N identical particles 
interacting by pair potentials V(rij) is given by 

JV N 

# = E\/^ +m2+ E ( L1 ) 

i— 1 

We use the notational simplification p= || p || , r = ||r||, or = || r j — r j || , whenever no 
ambiguity is introduced by so doing. Many approaches to such many-body problems 
for identical particles employ the very powerful constraint of permutation symmetry to 
generate their reduction to a two-body problem with a Hamiltonian H whose spectrum 
is used to approximate the many-body energy eigenvalues or to generate a lower 
energy bound. This reduction may be effected in various ways, which leads to the 
problem of finding the most effective reduced problem, the one which would provide 
the highest lower bound. In one analysis [6] involving pseudo-fermions (where the 
necessary permutation antisymmetry is carried entirely by the spatial part of the wave 
function), an optimization is considered over a real parameter which characterizes the 
degree of orthogonality of the matrix B that defines the relative coordinates. For boson 
systems, an orthogonal B is best possible; one such choice is the Jacobi coordinate 
system that we shall use in Section 2 below. 

For the boson problem, perhaps the most immediate reduction is what we have 
called the simple or N/2 bound based on the equality (H) = (H 2 ), where 



H - N 
H 2 -- 



\jv\ + m 2 + yjpl + m 2 + (N — l)V(r 12 ) 



(1.2) 



The N/2 bound is then the bottom E 2 of the spectrum of the scaled two-body 
Hamiltonian H 2 . We have explicitly for this bound 



E>E^ /2 = Nmi (j,, 



V> • (1-3) 



If this reasoning is applied to the Schrodinger harmonic-oscillator problem, one finds 
for large- N that E^, 2 — ► E/\/2, whereas a reduction based on Jacobi coordinates 
[11] yields E L — E. We note parenthetically that the N/2 bound is equivalent to 
using a non-orthogonal coordinate system consisting of a centre-of-mass coordinate 
and N — 1 pair distances [7] . Similarly, one obtains dramatic improvement over the 
N/2 lower bound if Jacobi coordinates are used for the Salpeter harmonic-oscillator 
problem [8] . We have obtained improved lower bounds for potentials which are convex 
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transformations V(r) = g(r 2 ) of the oscillator [9], and also, by very special reasoning, 
for the gravitational potential [10], V(r) — —v/r, v > 0. In the present paper we look 
for good lower bounds that are valid for general attractive potentials, for example, of 
the form V(r) = —v/r + br, v > 0, b > 0. 

Since the spectrum of the semirelativistic many-body Hamiltonian H can be 
characterized variationally, it is straightforward to find upper energy bounds with the 
aid of a suitable trial function. The principal difficulty is to find a good general lower 
bound. We achieve this for TV = 3, and for the case m = 0, N = 4. These partial results 
then allow the construction of corresponding lower bounds based on reductions of the 
many-body problem respectively to scaled N = 3 and N = 4 systems. A formulation 
that unifies these results and all the known earlier partial results may be expressed as 
a lower-bound conjecture, which then becomes a theorem for each case that is proved. 



2. Lower-bound conjecture 

We first consider a model iV-body Hamiltonian. This model has been constructed so 
that it essentially yields the corresponding nonrelativistic result in the limit m — > oo. 
We are guided in the first instance by the centre-of-mass identity and inequality [11] 

JV N / N \ 2 N 

i— 1 l—i<cj \ i / l—i<j 

For the corresponding semirelativistic problem, we lose this transparent algebraic 
inequality and must instead recover whatever can be proved to be true on the average. 
In a nutshell, this is the technical difficulty we must face in this paper. The Schrodinger 
iV-body Hamiltonian Hs with the centre-of-mass kinetic energy removed and % = 1 
is therefore given by 



i=»<j 



N 1 



2mN 



(Pi - VjY + v( rij ) 



(2.2) 



In Jacobi coordinates [p] = B[r], where B is an orthogonal N x N matrix with first 
row having entries all equal to l/\/N, p2 = (ri — r2)/\/2, and conjugate momenta 
[tt] = (B*) _1 [p] = B[p], the equality in (2.1a) may be re-written simply 



N N 



£p?=«? + E*?- (2-16) 



i=l i=2 



Meanwhile, if ^>(p2, P3, ■ ■ ■ , Pn) is a normalized translation-invariant iV-boson wave 
function, we have [10, Eqs. (6) and (7)]: 

(vl/,^ 2 *) = (*,7rl*) , (*,^*) = i = 2,3,.... (2.3) 

We note parenthetically, for future reference, that with Jacobi coordinates we have 
the following explicit expression for p^y: 



7Ti IN — 1 
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When either the kinetic energy is a quadratic expression, as for all Schrodingcr 
problems [11], or if the potential V(r) is the harmonic oscillator V(r) = kr 2 [9], 
then these relations play a key role in the construction of a lower-bound model. Our 
purpose here is to make a reduction for the Salpeter problem and general V(r), that 
is for problems for which neither the kinetic energy nor the potential energy has 
a simple quadratic form. We focus our attention on the kinetic energy since any 
progress made here would be potential independent. With these goals, the model 
iV-body Hamiltonian we have constructed is given by 



N 

E 

l=i<j 



'7" 1 (Pi - Pi) 2 + 



2m 

N -\ 



+ V(r l0 ) 



(2.5a) 



or, equivalcntly, 



N 



H c = £ 



l=i<3 



7(p* - p,) 2 + MV) 2 + v( rij ) 



(2.5b) 



where 7 = (^) = ^N(N — 1) is the binomial coefficient. In the Schrodinger limit 
m —* 00, we find H c — » mN + H$, where H$ is exactly the corresponding Schrodingcr 
A^-body Hamiltonian with the centre-of-mass kinetic energy removed, given in (2.2). 
Meanwhile, for the special case N = 2 of the semirelativistic problem itself we recover 
the well-known 2-body Salpeter Hamiltonian 



H 




+ m 2 + V(n 2 ). 



(2.6) 



If we use new conjugate coordinates, we may write r = 
II (Pi — P2)/2||, and then we have from (2.6) 

H = 2y / p 2 +m 2 + V(r). 



ri2 and p 



(2.7) 



The idea is eventually to obtain an TV-body lower bound which is the bottom of the 
spectrum of a scaled version of (2.6), namely 



H = /V V + m 2 + 7 V(r), /3,A, 7 >0. 



(2.8) 



Meanwhile, the Salpeter Hamiltonian H itself is given by (1.1). We now suppose that 
* is a translation-invariant normalized boson trial function. We consider expectations 
with respect to 4' and we first observe that the permutation symmetry of W implies 
the equality 

(H C ) = (H), where = N, A = >y=±N(N-l). (2.9) 

With these explicit values for the parameters {[3, A, 7} in H, wc are now able to 
formulate the central idea of this paper explicitly, namely 



Conjecture 



(H) > (H). 



(2.10) 
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This implies the following explicit conjectured lower energy bound 



E>E^ = iVinf V, 



N 



p 2 +m 2 ) + - V(r) 



V- • (2-11) 



We can recover all earlier sharp bounds from this expression. We immediately recover 
the Schrodinger bounds [11] in the to — > oo limit (2.5). If we now assume (2.11) is 
true as it stands for to > 0, and V(r) = vr 2 , we recover our earlier bounds [8] for 
the semirelativistic oscillator. For m > 0, and V(r) = —v/r, we recover our earlier 
sharp bounds for the gravitational problem [10]. Meanwhile, the bounds we prove in 
the present paper establish a wider range of validity for this conjecture. For example, 
our Theorem 3 below establishes (2.11) for to > and N = 3 in dimension d = 3; 
and Theorem 4 establishes the case m = 0, N = 4. At present we know of no counter 
example. 

If we compare (2.5b) with (1.1) we see that the expectation of the difference may 
be written 

(H - H c ) = {H -U) = (S(m,N)}, (2.12) 

where 



N i 9 N rw~\ 

S(m,N) = J2^+m 2 j^j ]T ^/ — ( P! - Pj ) 2 +ml (2.13) 

i— 1 l—i<j 

All our lower-bound results follow from the positivity (strictly speaking, non- 
negativity) of (6(m,N), when this can be established. We consider immediately the 
case {to = 0, N = 2}: the kind of reasoning we use in this case is generalized for 
the other cases. The approach we adopt is to think of the mean-value computation 
in momentum space where the momentum vectors Pi are multiplicative operators: 
these vectors form geometrical figures whose edges are the corresponding norms ||pi||; 
mean values (||pj||) = d arc considered last. For example, with N = 2, the three 
vectors {pi, P2, Pi — P2} form the sides of a triangle. The observation that, as a 
consequence of the triangle inequality and boson symmetry, the largest possible value 
for (||pi — p 2 ||) is 2d, then establishes positivity in this case. For to > the argument 
must be adjusted accordingly. We shall consider this point in more detail in Section 5 
below, for the more interesting case N — 3 and to > 0. In order to prepare for what 
might be called 'stochastic geometry', we consider first N = 3 and to = 0, although 
this is a special case of the more general problem to > discussed later. As we have 
remarked above, for the corresponding Schrodinger problem for general V(r), or for the 
Salpeter harmonic-oscillator problem with V(r) = kr 2 , a quadratic form is involved 
cither in the kinetic- or the potential-energy term: for both of these problems, the 
conjecture follows as a result of the general quadratic mean-value identities (2.3) in 
Jacobi coordinates. For the Salpeter problems with general V, which is the subject 
of the present paper, the quadratic expressions (in momentum space) always appear 
inside the square-root sign, so these identities do not immediately apply. The general 
inequality (|jp||) < (||p|| 2 ) 5 docs not remove this difficulty 
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3. Proof in the case m = 0, N = 3 



We have the following definition from (2.13): 

5(0,3) = HpxII + ||p 2 || + HpsII - -i= (||pi - p 2 || + ||pi - Pall + ||p 2 - Pall) • (3-1) 

(S(0, 3)) = ^llpiH + ||p 2 || + HpsII - ^= (||Pi - P2 || + ||pi - Pall + ||P2 - Pall)) • 

(3.2) 

We note that 5(0, 3) itself is negative for the choice p 2 = Pi ^ and P3 = 0. 
However, this does not happen on the average. We have: 



Theorem 1 (5(0,3)) > 0. 



Proof We know by boson symmetry that 

(l|pi||> = (l|P2||) = (||Pa||):=fc (3.3) 

and 

(||P1-P2||) = (||P1-P3||) = (||P2-P2||):=9. (3.4) 

We may think of the {p^}, and their differences, as vectors, since they are used in 
momentum space where they become multiplicative operators. The six vectors in 
(3.1) are the six edges of a pyramid in 3? 3 ; the norms, ||p,|| and \\pi — pj\\, are 
the corresponding lengths of these six pyramid edges. The permutation symmetry 
of the wave function implies the equalities (3.3) and (3.4). The mean difference 
(5(0,3)) is clearly smallest when the origin of the vectors {pj} is at the centroid 
of the triangle formed by the differences {p.; — pj}. In this minimal case we know 
from elementary geometry that q = \^3k; consequently, (5(0,3)) = 0. It follows that 
in general (5(0, 3)) > 0. This completes the proof for the case m — 0, N = 3. 

□ 



4. Proof for the case m = 0, N = 4. 



We consider the case N = 4 and m = in (2.13). The six differences {pi — Pj} 
form a tetrahedron. The average lengths q — (||pj — Pj||) are equal and force the 
tetrahedron to be regular. Meanwhile, the four mean lengths k = (||Pi||) are again 
equal. This symmetry occurs when the p-origin is at the centroid of the tetrahedron, 
of, say, height h. For such a tetrahedron we have 



h 



2 /3 
—q and k = \ -q. 

3 V 8 



(4.1) 



We may therefore write 



( M0.4))=4(|| Pl .||; - (,(l)Jl { \\ P] P2 | 



4fc - i\ -q = 0. 



(4.2) 



Thus we have: 



Theorem 2 (5(0,4)) > 0. 



Binding energy of semirelativistic N -boson systems 
5. Proof in the case m > 0, N = 3 
We consider 



6(m,3) = (|| Pl || 2 + m 2 )^ + (||p 2 || 2 + m 2 )* + (||p 3 || 2 + m 2 ) * 

- (illPi - P2|| 2 - (i|| Pl - p 3 || 2 + m 2 )^ - (i|| P2 - P3 ||2 + m ?)h (5.1) 

and 

(5(m,3)) - ((|| Pl || 2 + m 2 ) h + (||p 2 || 2 + m 2 ) ' + (||p 3 || 2 + m 2 ) * 

- (|||pi - P2 || 2 +m 2 )^ - UHpi - P3 || 2 + m 2 )^ - (i||p 2 -p 3 || 2 + m 2 )^. (5.2) 



Theorem 3 (6(m,3)) > 0. 

Proof The three vectors Pi, i = 1,2,3, and their three differences p^ — Pj form 
six edges of a pyramid in 3? 3 ; the norms, ||pj|| and ||pi — Pj||, are the corresponding 
lengths of these six pyramid edges. We now denote by T the triangle formed by the 
three difference edges {||Pi — Pj||}- For convenience, we shall think of T as lying in 
a horizontal plane and denote by P the top vertex of the pyramid; without loss of 
generality, we shall speak of P being above T. We let C be the point in the plane 
of T vertically under P. We now pick the vertex of T which contains pi, and call 
this V\. In the plane of T we construct a line from V\ that is perpendicular to CV\ 
and of length m, ending in the point U\. We then join U\ to P and observe that 
PV\U\ = 7r/2. Similar constructions arc now made with the other two vertices V 2 
and V 3 of T; the three line segments XJjVi are chosen to 'flow' in the same way round 
the pyramid axis CP. In fact, a new pyramid is formed by the three lines PUi. By 
permutation symmetry we have that (\PUi\) = k and (|Cf7j|) = q where i = 1,2,3, 
and moreover 



((|| Pj || 2 + m 2 )^ :=k, z = l,2,3, (5.3) 



and 




p,|| 2 + m 2 ):=q, i,j = 1,2,3, i^j. (5.4) 



Since the position of P which minimizes k is C, and symmetry is obtained on the 
average, we conclude by elementary geometry that k > q. This inequality completes 
the proof of Theorem 3. 

□ 



6. Application to N > 3 

For N > 3 we can deduce a stronger lower bound than that provided by the N/2 
bound, based on the result of Section 5. If E and ^ are the exact energy and 
corresponding wave function, we have E = (^f, ff VP) and therefore, by boson symmetry 
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and Theorem 3, we have 
TV 



E= — (tf , [{p\ + m 2 )3 + {pi + m 2 )i + (pj + m 2 )i 
+ ^Y^(V(r 12 ) + V(r 13 ) + V(r 23 )) * 

/ r 1 . . , n - 1 i 

> w I *, 

> tv [ 



/ 1 9 ?\i JV - 1.,, \ 

(gPi2 + "» ) 2 +— 2~ V(n 2 ; 



(-p 2 + m 2 ) 2 +__y( r ) 



where r = ri — r 2 and p = |(pi — p 2 ) = pi 2 - Thus we have, for TV > 3, m > 0, and 
||V(r)||=l: 



Theorem 4 



E > E 



N/3 



TVinf V, 



2 2 

-p + m 



TV — 1 



V(r) 



1>\. 



(6.1) 



In similar fashion, we can relate the TV-body problem for TV > 4 and m = to a 
reduced 4-body problem based on Theorem 2. Specifically, we have for iV > 4, m = 0, 
and ||V>(r)|| = 1: 



Theorem 5 



E>Ek /4 = NM (V, 



;|p|| + ^rr^^W 



(6.2) 



Theorems 4 and 5 summarize the principal results of this paper. 



7. The linear potential V(r) = r with m = 

The lower bounds we have found all presume that the bottom of the spectrum of a 
scaled one-body problem can be found. For Salpeter Hamiltonians, this task itself may 
not be trivially easy, although more tractable than for the many-body problem. For 
the operator H = ||p|| + r in three dimensions, we have at our disposal the accurate 
value e = 2.2322, for example, from the work of Boukraa and Basdevant [12] (the 
linear potential has also been considered by Pirner and Wachs [13] in an application 
to quark systems). By elementary scaling arguments we therefore have for the one- 
body problem 



H = ap + br 



E(a,b) = (ab)?E(l,l) = (ab)*e, a,b>0, e = 2.2322. (7.1) 



This relation will generate all the energy lower bounds for TV-body problems with this 
potential. We shall use the notation £^y 2 , ^/v/3> an( ^ ^v/4> f° r tnc l° wcr bounds 
given by equations (1.3), (6.1), and (6.2), and E c for the conjectured bound (2.11). 
The formula (7.1) then allows us to derive formulas for these energies. Explicitly we 
find: 

1 

/ iV — 1\ 2 

E N /2 = N — — e, TV > 2 (7.2a) 
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/ AT — 1 \ 2 

E L N/3 = N e, N > 3 (7.26) 

E n/4 = N ( 3iN ~ 1)2 ) 4 e, JV > 4 (7.2c) 



^ = TV I (jV oAr 1)3 I e, N > 2. (7.2d) 



27V 

In order to find an upper bound, we follow Ref. [9] and use a Gaussian wave 
function, which wc write initially in the form 

$(p 2 ,p 3 ,..., PA r) = Ccxp --^^ 2 U[]^), C= — , (7.3) 

V t=2 / i=2 VTl" 4 / 

where the constant C is chosen to ensure the normalization of each radial factor 
4> on L 2 ([0, oo), r 2 dr). The boson symmetry of the trial function allows us to write 
E < E^ = ($, $) , where we have 

£* = ($, [TVHpjvH + 7 ^(||n - r 2 ||)] <f>) . (7.4) 

The identity (2.4) and the lemma proved in [8] (which allows us to remove the operator 
term m) imply 

E u g = N^j Uttjv II + ~fV(V2p 2 )$j . (7.5) 

The permutation symmetry of the Gaussian function in the relative coordinates and 
the factoring property allow us to replace ttjv by 7r 2 = V2p. We write the conjugate 
variable to p as r e V%P2, so that V(r) = r, and the wave function becomes 4>{r). By 
introducing an additional scale parameter a > 0, we then hnd 




Since the Gaussian radial function <f>(r) is form invariant under the 3-dimensional 
Fourier transformation, we have the equality 

(P> = (r> = -^=. 

Since the minimum of the form a/a + [3a over the scale a > is 2(a[3)i , we arrive at 
the following explicit formula for the Gaussian upper bound: 

= 4N 4 ' N * 2 - (7 ' 7) 

We can immediately test this formula for the case N = 2 to obtain E^ = 3.19154, 
which is to be compared with the accurate numerical value E = \[2e = 3.1568. More 
generally, we exhibit in Table 1 ratios R x = E v g /E%, where X is N/2, N/3, N/4 or, for 
the conjectured lower bound, R c = Eg /E^. The percentage error in the determination 
of the energy by the bounds is approximately 50 (R — 1)%. The monotonic behaviour 
of R with N follows from the 'distance' of N from the size of the sub-system whose 
lower bound is best possible; if the conjecture were true, the quality of the lower bound 
would be the same for all N. 
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Table 1. Ratios of upper to lower energy bounds Rx = E]j /E x , where 
X = N/2, N/3, N/4; the ratio for the conjectured lower bound is R c = E^/E^. 





TV = 2 


TV = 3 


TV = 4 


TV = 5 


TV = 6 


TV = 10 


TV —> oo 


Rn/2 


1.011 


1.08639 


1.11886 


1.13706 


1.14872 


1.17104 


1.20229 


Rn/3 




1.011 


1.04121 


1.05815 


1.069 


1.08977 


1.11886 


Rn/4 






1.011 


1.02745 


1.03799 


1.05815 


1.08639 


Rc 


1.011 


1.011 


1.011 


1.011 


1.011 


1.011 


1.011 



8. Conclusion 

If a system of TV identical particles is bound together by attractive pair potentials, the 
Hamiltonian H has TV kinetic-energy terms and 7 = (^) potential terms. If the kinetic 
energy of the centre-of-mass can be subtracted off, then the number of kinetic-energy 
terms is reduced by one, and we would expect to obtain an expression of the form 
E = {H) = ((TV — 1)K + jV). The TV-body energy E is then bounded below by the 
lowest energy £ of a 'reduced' one-body operator of the form H = (TV — 1)K + -fV; 
if the boson-symmetry requirement of the TV-body wave function is not too stringent, 
then this lower bound is at the same time a good approximation. This story is realized 
exactly for the nonrelativistic problem [11]: for the special case of the harmonic 
oscillator, £ yields the exact energy E of the many-body system. The reduction details 
depend on the quadratic form of the nonrelativistic many-body kinetic-energy operator 
and the identities (2.3) for quadratic expressions in Jacobi relative coordinates. 

For the semirelativistic counterpart, one generally loses the quadratic form in 
H and, along with it, the immediate reduction. An alternative reduction to the 
Hn/2 Hamiltonian is always possible and is important theoretically, but the resultant 
lower energy bound is weak. A quadratic form is returned to the potential in H in 
the special case of the harmonic oscillator, and this yields [9] a very sharp bound 
on the energy, though not now the exact solution, except in the Schrodinger limit 
m — > 00. For general pair potentials, we have constructed a new Hamiltonian H c that 
is translation invariant, both in coordinate and momentum space, and which reduces 
to the usual two-body Hamiltonian for TV = 2, and generally to Nm + Hs in the large- 
m limit, where Hs is the Schrodinger Hamiltonian with the centre-of-mass kinetic 
energy removed. A reduction (H c ) — (TL) > £ of H c to a one-body Hamiltonian Tt 
immediately follows. This is useful for the study of the many-body Hamiltonian H 
whenever it can also be established that (H) > (H c ). We conjecture that this is always 
true. In the present paper we have proved the conjecture for TV = 3, and for N = 4 if 
m = 0; it is also true for the harmonic oscillator, and generally in the largc-m limit. 
For the case of static gravity V(r) = — 1/r, the conjecture yields the identical result 
to the energy bound we have established by a completely different argument, valid 
specially for this potential [10]. 
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